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ABSTRACT 

We present a morphological analysis of the feature-rich 2MASS LMC color- 
magnitude diagram, identifying Galactic and LMC populations and estimating the 
density of LMC populations alone. We also present the projected spatial distributions 
of various stellar populations. The conditions that prevailed when 2MASS observed 
the LMC provided lOo limiting sensitivity for 7 < 16.3, < 15.3, i^i- < 14.7. 
Major populations are identified based on matching morphological features of the 
color-magnitude diagram with expected positions of known populations, isochrone 
fits, and analysis of the projected spatial distributions. 

2MASS has detected a significant population of asymptotic giant branch (AGB) 
stars ( > 10^ sources) and obscured AGB stars ( > 2 x 10-^ sources). The LMC 
populations along the first-ascent red giant branch (RGB) and AGB are quantified. 
Comparison of the giant luminosity functions in the bar and the outer regions of the 
LMC shows that both luminosity functions appear consistent with each other. The 
luminosity function in central (bar) field has well-defined drop-off near Ks = 12.3™ 
corresponding to the location of RGB tip; the same feature is seen in the luminosity 
function of the entire LMC field. 

Isochrone fits for the corresponding giant branches reveal no significant 
differences in metallicities and ages between central and outer regions of the 
LMC. This may be evidence for strong dynamical evolution in the last several 
gigayears. In particular, the observed LMC giant branch is well-fit by published 
tracks in the CIT/CTIO system with a distance modulus of/j = 18.5±0.1, reddening 
Eb-v = 0.15-0.20, metallicity Z = 0.004+5 [JjJ? and age 3-13 Gyr. Analysis of 
deep 2MASS engineering data with six times the standard exposure produces similar 
estimates. 



Subject headings: astronomical data bases: surveys — galaxies: luminosity function, 
mass function — Magellanic Clouds — galaxies: stellar content — infrared: stars 
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Introduction 



Large and homogeneous data sets of near-infrared photometry for the entire LMC, a 
by-product of large-scale infrared sky surveys such as 2MASS ( [Skrutskie a/. 1997| ) and DENIS 
( |Hpchtein et al. 19^7] ), have become available to astronomical community only recently. Cioni et 



al. (2000) introduced the DENIS Point Source Catalog towards the Magellanic Clouds. Here, we 
present the LMC data from The Two Micron All Sky Survey (2MASS). 

The 2MASS has observed the entirety of the Large Magellanic Cloud and much of these 
data are included in the recent second incremental data release. Empirically, the photometry has 
signal-to-noise (SNR) ratio 10 at J, H, magnitudes of 16.3, 15.3, 14.7, respectively, slightly 
better than the nominal survey limit. At these limits, we can observe all of the thermally-pulsing 
asymptotic giant branch (AGB) populations and part way down the red giant branch (RGB). The 
red clump, representing helium burning giants, lies ~ 2 mag below the sensitivity limit of these 
data. The extinction in near-infrared (NIR) is small throughout the LMC and negligible on average 
everywhere but the inner degree of arc. Together with the high quality of 2MASS photometry 
(Oot ~ 0.03"), overall zero-point stability (better than 0.01'") and with reliable identification of 
LMC populations, the survey is ideal for studies of spatial structure of the LMC or its evolution 
(see the companion paper, Weinberg & Nikolaev 2000; hereafter WN). 

We describe data selection, cross -correlation with a few well-known populations, and 
comparison of the 2MASS color system in §^ We present a morphological analysis the 
color-magnitude diagram in §|3|. In particular, we make identifications of the Galactic and LMC 
populations corresponding to all features in the color-magnitude diagram, correlate these with 
spatial distribution, and estimate the density of LMC populations alone. The LMC giant branch 
is well determined and we separately identify the AGB, first-ascent red giant branch tip (TRGB), 
and the carbon star sequence. The luminosity function of RGB and AGB populations is derived 
(§Q) and compared with the galactic giant-branch luminosity function. We explore the feasibility 
of determining the spatial dependence of metallicity using giant-branch morphology. Finally, §^ 
summarizes our results and discusses the implications and opportunity for future study. 



2. Data 

Our LMC field is ~ 250 sq. degrees and covers the range from 4''00'" to 6^56'" in right 
ascension and from —77° to —61° in declination (coordinates in J2000.0). The initial sample of 
7,092, 894 sources is drawn from the Working Survey Data Base and includes possible artifacts, 
such as filter glints and diffraction spikes from nearby bright stars, source confusion, and detection 
upper-limits. The known contaminants and flux statistics are well-characterized and identified 
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during processing ( JCutri et al. 1999| ). Eliminating artifacts and requiring detections at all bands 
with Om < 0.11™ (SNR > 10) leaves 1,246,304 stars. Unlike the released catalog, these data 
contain multiple apparitions sources because of scan overlaps. We identify the multiple entries 
based on 1) spatial proximity (|Ar| < 2"), and 2) matching -photometry (|Ar| < 5o). Note that 
this procedure differs from that of the 2MASS catalog release (see Cutri et al. 1999). Our final 
sample contains 823,037 sources. 

The spatial distribution of these stars is shown in Figure [I]. The figure shows major structural 
components of the LMC, the bar and disk, immersed in the field of Galactic foreground. The 
gradient of the foreground sources (the direction to the Galactic center is indicated by an arrow) 
distorts the isopleths of the outer LMC disk. The source density near the optical center of the bar 
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Fig. 1. — Distribution of 2MASS sources in the LMC field. Contour levels are labeled in units of 
10^ deg^^. Arrow points in the direction to Galactic center. 

(at a = 5''24'", Sy2000 = — 69°44') exceeds 3.6 x 10^ deg ^. The expected mean separation among 
sources at this density, ^ = 1/ y/p ~ 20", is much greater than 2MASS resolution and therefore 



-4- 



confusion should not be significant (see also Wood 1994). As a separate check, we have examined 
the source counts as the function of magnitude in dense central fields and in more sparse fields in 
the outer LMC. The distribution of source counts as a function of magnitude have similar shapes 
in both dense and sparse regions, consistent with low confusion. 

We have cross -correlated our sources with the database of long-period variables (LPV) from 
Hughes & Wood (1990). Their sample includes 376 large-amplitude Mira-like variables and 224 
smaller amlitude semi-regular variables. We find 370 (98.4%) and 224 2MASS counterparts, 
respectivel)!^. Three of the 'missing' large- amplitude variables are present in the raw 2MASS 
data, but are degraded by artifacts (two diffraction spikes and one blend). Of the remaining three, 
two are matched by stars of the appropriate magnitudes within 1" radius, and only one does not 
have any match to within 20" radius. All 134 Wolf-Rayet stars in the LMC ( |Breysacher 1999[ ) 
have been observed by 2MASS ( |van Dyk g?^ a/. 1999| ). 



2.1. Ks vs. K Photometry 



The 2MASS photometric system is similar to CIT/CTIO system ( plias et al. 1982D , except 



that it uses Ks band (2.00 — 2.32;um) rather than K band. The Ks ('K-short') bandpass is described 
by Persson et al. (1998). It was designed to reduce the ground-based thermal background. The 
transmissivity curve for the filter is given in Persson et al., who also compared Kqit with Ks 
photometry for a set of solar- type stars and red standards (see their Tables 2 and 3, respectively). 
Based on their data, the difference K — Ks shows no significant systematic trend in the color range 
< 7 — ^ < 3. The strongest trends follow from the presence of CO-band absorption, which 
affects the K filter more than the Ks filter. The absolute value of the difference |^ — is less than 
0.05'" and we will assume Ks = Kcit in comparing CIT/CTIO system-based stellar sequences 
with 2MASS data. 



2.2. Interstellar Reddening 

One of the advantages of 2MASS as compared to optical surveys is low interstellar reddening, 
since extinction at 2;um is approximately 10 times smaller than in V. The values of interstellar 
reddening Eb-v found in the literature^ fall in the range between 0.08 ( fVIateo & Hodge 1987D 
and 0.20 ( [Harris et al. 1997| ). The distribution of reddening from Harris et al. (1997) has 



The search radius extends to 5" from the listed positions of sources, 
^foreground plus LMC internal reddening 
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non-Gaussian tail to high values. Greve et al. (1990) have reported values as high as Eg-v = 1- 
found from their investigation of dust in emission nebulae in the LMC. Bessell (1991) has 
summarized reddening determinations from photometry, stellar polarization and HI column 
densities to derive foreground and intrinsic mean reddening in the Clouds. He obtained typical 
LMC internal reddening of 0.06 (with substantial variations), and foreground reddening in the 
range 0.04 — 0.09. Galactic foreground reddening can be surprisingly large in the outer regions of 
the LMC: Walker (1990) reported Eb-v = 0.18 ±0.02 at NGC 1841, about 15° from the optical 
center of the LMC. On the other hand, in the cluster GLC 0435-59 (Reticulum), 11° from the 
LMC center, the reddening is only Eb-v = 0.03 ( [Walker 1992| ). Zaritsky (1999) has indicated 



that reddening for F and G stars in the LMC is much less than that for OB stars and derived the 
average < Eb-v >= 0-03 for late-type stars in the disk. 

In the present study, the data are not dereddened. Rather, each diagram shows the direction 
and magnitude of the reddening vector for a specified value of Eb-v- The reddening vector is 
based on relations from Koomneef (1982): 

AK = 0.mEB-v; Ej-k = 0.651Eb-v, 

fori? = Av /Eb-v = 3.1. Information about the reddening may be obtained directly from 2MASS 
data from the analysis of the color-color diagram (see below). 



3. Analysis of the Color- Magnitude Diagram 

Figure shows the color-color {J-H vs. H- K^) diagram of the LMC for 823, 037 2MASS 
sources selected in §|[ The diagram shows relatively few distinct features. The most prominent 
among them the extended 'arm' of the thermally-pulsing AGE stars (TP-AGB) in the upper 
right corner. Typical colors of LMC M giants are in the range 0.2 < H — ^ 0.3, J — H > 0.8 
( progel & Blanco 1990| ). Most stars on the extended arm, with colors redder than J — Ks = 1.6 are 
carbon-rich stars ( Pughes & Wood 1990D . Of those with the extremely red colors, J — Ks > 2.0, 



many probably possess dusty circumstellar envelopes. The sample contains ~ 2000 such sources. 
Their locus is consistent with the track following the reddening vector for H — Ks> 1 .0"\ Figure ^ 
also shows fiducial color tracks for both giants and dwarfs from Wainscoat et al. (1992; hereafter 
W92). The two sequences overlap near H — = 0.15, J — H = 0.5, since NIR colors of late G — 
early M type dwarfs are the same as colors of late F — early K type giants. 

Because of the small number of features and the general compactness of the color-color 
diagram, its usefulness in discriminating the major populations is limited, especially in the overlap 
region, 0.5 < J — < 0.8. However, some LMC populations which occupy distinct regions in 
the diagram can still be identified based of their location in the color-color plot (see Figure 0). 
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Fig. 2. — Left panel: Color-color diagram of the LMC field. Contour levels are logarithmic, from 2 
to 6.5, spaced by 0.5. Diagonal lines are lines of constant J — Ks values (marked). Color sequences 
of dwarfs (solid line) and giants (dashed line) from W92 are shown. Reddening vector for 
Eb-v = 1-0 (indicated by arrow) is based on relations from Koornneef (1982), assuming i? = 3.1. 
Right panel: Color-color diagram showing approximate positions of some LMC populations. 
Shaded area corresponds to Wolf-Rayet stars (Breysacher 1999); region outlined by dashed lines 
encompasses known LMC B[e] stars (Gummersbach et al. 1995), open squares show individual 
observations of four LMC protostars (from Westerlund 1997). W92 fiducial colors of dwarfs (thick 
solid line) and giants (thick dashed line) are indicated. Region occupied by carbon stars in the 
sample of Costa & Frogel (1996) is shown with solid lines. Dotted line corresponds to 7 — ^5 = 1.6. 
Reddening vector is drawn for Eg-v = 1 -0. 

In particular, the color-color diagram is quite useful in identifying candidate sources with large 
infrared excess, such as young protostars, cocoon stars, or obscured AGE carbon stars. 

The color-color diagram may be used to determine the reddening distribution. The giant 
population forms a tight branch in near- infrared colors. The reddening direction nearly coincides 
with J — Ks and therefore the distribution in 7 — ^ of a sample from narrow color interval in 
J — H along the giant branch provides a sensitive diagnostic for reddening by dust. We considered 
a sample of sources in the range 0.78 < J — H < 0.85. The resulting shift of the peak is 
A(7 — ^5) < 0.03 (A{J — Ks) < 0.06) outside (inside) of the central region, suggesting only minor 
reddening on average on scales larger than 0.1 square degrees. 
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Fig. 3. — Left panel: Color-magnitude diagram of the LMC field. The density levels are 
logarithmic, from 2 to 6, spaced by 0.5. The reddening vector corresponds to Eg-v = 1.0. Right 
panel: The same diagram with highlighted 12 regions discussed in text. The regions correspond to 
major features of the CMD. 

The color-magnitude diagram (CMD) presented in Figure |] reveals a wide variety of details. 
Our goal is reliable identification criteria for LMC stellar populations based on their positions in 
the diagram. The CMD is hand-shaped with vertically-stretched 'fingers' (e.g., at 7 — Ks colors of 
0.4, 0.6, 1.1) due to varying distance modulus for both Galactic and Magellanic sources. We have 
identified 12 regions shown in Figure ^ that highlight distinct features of the CMD. The regions 
are marked A through L and enclose 99.7% of the 823, 037 sources in the field. To identify stellar 
populations in each region, we use a combination of several techniques. The Galactic foreground 
contribution is modeled by a synthetic CMD based on the tabulated near-infrared model of W92. 
The LMC populations are identified based on the infrared photometry of known populations 
found in the literature. In addition, we do a preliminary isochrone analysis, where we match the 
features of the CMD with Girardi et al. (2000) isochrones to derive the ages of populations and 
draw evolutionary connections among the CMD regions. The details of our population matching 
procedure are given in § p.l| . In addition, we use the spatial density distributions of sources in each 
region to better discriminate local and LMC populations. The spatial distributions for each region 
are shown in Figure ^. In each panel of the figure, we plot source density contour levels of 15, 30, 
60, 120, 240, 350, 480, 960, 1920, 3840, 7680, and 15000 deg"^ j^e contour levels are selected 
both to show the underlying density profiles with maximum details and facilitate the comparison 
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Fig. 4. — Spatial density distributions of sources in CMD regions. Letters correspond to the 
regions introduced in Figure ^. The distributions are kernel smoothed with a Gaussian kernel. The 
same sequence of contour levels is used in each panel (see text). Contour level of 120 deg^^ is 
highlighted. 
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among different population densities. However, due to strong variations in relative density in the 
CMD, not all contour levels in this sequence are displayed: in panels E and I the lowest density 
contour corresponds to 60 deg^, in panels B and C the contours start from 240 deg^, and in very 
dense Region D the lowest contour level is 960 deg^^. 



3.1. Identifying Stellar Populations 

The initial analysis of the CMD regions has two parts: 1) use of the spatial density distribution 
to estimate the location (Galaxy or LMC); and 2) use of the theoretical colors/isochrones to derive 
the properties of the population, such as age, approximate spectral class and distance modulus. 
Here we describe the procedure of identifying the populations, before examining each region of 
CMD in detail. 

The CMD of the LMC field (Fig. ||) contains both Galactic and LMC populations. To 
quantify Galactic foreground, we use near-infrared model of W92, based on 8-25 micron point 
source counts. Galactic model of W92 has five structural components: exponential disk, bulge, 
stellar halo, spiral arms and molecular ring. The main contribution to the Galactic source density 
toward the LMC (/ = 280.5; b = —32.9°) is the exponential disk. The source density due to 
other Galactic structural components combined does not exceed 0.0005% in any region of the 
CMD. The luminosity function in W92 model is represented by a sum of stellar classes, allowing 
independent estimate of the contribution of each class to the CMD. Each class of source is 
assumed to have a Gaussian distribution. 



N{M) oc exp 



2o2 



To model the Galaxy, we use the first 33 classes from Table 2 of W92 (Galactic dwarfs, giants 
and supergiants). The remainder (AGBs, planetary nebulae, etc.) are expected to give only a 
small contribution to source density and thus do not affect the CMD. In the Galactic model, we 
use reddening parameters from Rieke & Lebofsky (1985). Dust is assumed to follow a double 
exponential distribution, with the radial scale length of the disk and the scale height of 100 pc. 
We reduced the magnitude dispersion a of each stellar class by ten to accurately represent the 
CMD. Reducing a results in a spiky differential luminosity function but this does not affect 
our application and the cumulative luminosity function remains well-approximated. Given the 
granularity of the model, the agreement between original W92 luminosity function and ours is 
acceptable. Our synthetic 'foreground' CMD is shown in Figure |], along with the observed 
CMD of a Galactic field 0. The agreement between the model and observed CMD is good, with 



^To define our Galactic field, we combined three small fields, each w 0.6 sq. deg. These fields were selected at the 
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Fig. 5. — Comparison of (a) observed CMD of a Galactic field and (b) the synthetic Galactic CMD 
from W92 model, as explained in text. Densities in both diagrams are normalized to unity. The 
contour levels are logarithmic, from —0.5 to —3.5, spaced by 0.5. Similar plots for each stellar 
class and for each structural component of the Milky Way allows unambiguous determination of 
populations responsible for observed CMD features. 

Conclusions about the LMC populations in the CMD are made based on isochrone fitting 
or on empirical matching of populations found in the literature to features of the CMD. We 
use theoretical isochrones from Girardi et al. (2000). These isochrones supersede Bertelli's set 



( iBertelli e?^ g/. 1994[ ), and use updated opacities and equations of state. The isochrones follow the 
evolution of low- and intermediate-mass stars (0.1 5M0 < M < IMq) from the main sequence 
up to the tip of the RGB or the start of the thermally -pulsing AGE. Distinct LMC populations 
are identified by matching morphological features of the CMD with colors of known populations 
from the literature. In particular, we use Cepheid colors from Madore & Freedman (1991), early 
M supergiant color sequence from Elias et al. (1985), and data on long-period variables from 



boundary of our LMC field where contamination from LMC source is minimal. At the time, other contiguous fields 
were not available. 
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Hughes & Wood (1990). This matching is purely qualitative and is only used as a supplement. 
The fiducial colors of Galactic giants and supergiants from W92 model are unsuitable for LMC, 
since LMC has a lower metallicity compared to the Milky Way. Moved to the LMC distance, 
/J = 18.5, the W92 giant branch provides a poor fit to the observed RGB (see Figure [7| below). 



3.2. Region A: Blue Supergiants, O Dwarfs 

These blue-colored sources are readily identified as early type Population I stars in the 
LMC. This group of stars is the evidence of recent (< 30 Myr) star formation. Plotting the 
theoretical evolutionary tracks in the CMD (Figure |^) confirms that the region is populated by 
blue supergiants and brightest dwarfs (ZAMS). Only the hottest and most massive dwarfs of types 
03 — 06 can be seen in the LMC at /j = 18.5. All other main sequence (MS) populations are 
too faint and fall below the limit at SNR = 10 imposed for this work. The supergiant population 
in the region are core helium-burning stars with masses 4 < M < 9Mq. These stars spend most 
of their post-MS lifetimes as blue or red supergiants QMaeder & Meynet 1989| ) looping between 



Regions A and H of the CMD (cf. Figure ^b). Region A encompasses stars which are at the blue 
tips of their blue loops. While crossing Regions B and C, these stars enter the instability strip and 



become Cepheids (see §§ |3]J, [J4| ) 



The spatial distribution of these objects also clearly indicates an LMC population. The 
distribution is rather clumpy with several richest OB associations outlining the location of spiral 
arms and brightest and largest HII regions (e.g., 30 Dor). The density concentrations in Fig. ^ 
are cons istent with the well-known superassociations and Shapley's Constellations (^M artin et 



al. 1976; [van den Bergh 198 1| ). These youngest populations do not trace the bar of the LMC, 
in agreement with de Vaucouleurs & Freeman (1973). Quantitative analysis of the distribution 
( [Weinberg & Nikolaev 20"00| ) puts the centroid of the population at a = 5''23'", 5 = -68°48', 
about 1° north of the optical center of the bar. 

The Galactic population of early dwarfs is readily seen in the CMD directly above Region A, 
blueward of Region B. The apparent magnitudes of these stars suggest a distance modulus between 
5 and 10 (r ~ 0. 1 — 1 kpc). In addition, this area of the CMD may also contain contribution from 
field blue stragglers and blue horizontal branch stars. 



3.3. Region B: Galactic Disk F — K Dwarfs, LMC Supergiants 

Region B is a vertically stretched band in the CMD with J — Ks = 0.2 — 0.5. This color cut 
isolates the main sequence tumoff of the halo (7 — ~ 0.3) and the disk (7 — ~ 0.4). The 
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spatial density distribution increases toward NE comer of the field (Fig. i.e. toward the 
Galactic center, and indicates a predominantly Galactic population. The vertical extent in the 
CMD indicates a wide range of distance moduli for these stars. Based on relative population 
abundances in our synthetic Galactic CMD, we conclude that these sources are disk dwarfs 
of spectral classes in the range from late F to early K. These stars account for ~ 90% of the 
foreground source density in the region. Their position in the CMD (see Figure ^ja) suggests that 
the dwarfs have distance moduli /j = 3 — 10 (r ~ 0.04 — 1.0 kpc). Galactic giants in the region are 
of types F — G, but their contribution to foreground source density is insignificant, smaller than 
5%. 

The distorted shape of the central isopleth in Figure ^ suggests presence of the LMC 
population in this region. The isopleth outlines the structure similar to the one seen in the 
central regions in Figure Note the overdensity near a = 5''40'", 5 = —69° and a = S'^BS'", 
6 = — 67°30', marking positions of superassociations IV and V, respectively ( [Martin et al. 197^ ). 



Based on the colors, the LMC component is comprised of young blue and yellow supergiants, 
corresponding to spectral types A — G. This population includes luminous blue variables and 
short period Cepheids {P < 50'^). Figure 0b shows the Cepheid sequence based on PL relations 
for LMC Cepheids ( [Madore & Freedman 1991D . Figure ^ also shows the colors of supergiants 



from W92 (Table 2). In addition. Region B contains the majority ( > 80%) of the known LMC 
Wolf-Rayet stars. Most LMC Wolf-Rayet stars have infrared colors in the range < 7 — < 0.5. 
Their numbers, however, are not significant to produce an observable effect in the CMD density. 



3.4. Region C: Disk K Dwarfs and K Giants, Young Supergiants in the LMC Bar 

Similar to B, Region C is stretched along the magnitude axis, indicating that the CMD feature 
is formed by sources at a range of distances. The colors of this population are in a tight range, 
A{J — Ks) ~ 0.3. Our synthetic Galactic CMD suggests that most (~ 70%) of the observed density 
in this region is produced by disk K dwarfs at < 9 (r < 600 pc). Disk K giants are also present 
in this region. Most of them have ~ 6 — 13 (r ~ 0.2 — 4 kpc). They contribute ~ 20% of the 
foreground density. The inspection of isochrones in Figure ^ suggests that Galactic giants in this 
region are in the evolutionary phase of red clump/horizontal branch stars. The intrinsic brightness 
and color of these stars {Mk = — 1.4±0.1, J — Ks = 0.6±0.1) make them natural candidates 
in this region. Because of the narrow magnitude range of the clump, the source distribution in 
Region C along the magnitude axis could help constrain the structure of the Galactic disk. 

The LMC population in Region C (seen in Figure ^fZ) is slightly older than youngest 
supergiants in Regions A, B. The central isopleths of the figure outline the bar of the LMC and 
show no overdensity at the positions of superassociations, seen in previous two regions. Most of 
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Fig. 6. — Part of the CMD showing Galactic and young LMC populations, (a) Galactic dwarf 
populations in 2MASS CMD. Solid line: theoretical isochrones for x = 7 Gyr, Z = 0.019, /j = 9.0, 
Eb-v ==0.1, representing intermediate/old Galactic disk population; dashed line: isochrones for 
X = 14 Gyr, Z = 0.0004, ^ = 9.0, Eb^v = 0.1. Fiducial unreddened dwarf colors from W92 (for 
/J = 9.0) are indicated by filled squares, with the approximate spectral types marked in the top 
axis. The reddening vector corresponds to Eb-v = 0-2. (b) Young LMC populations. Solid line 
shows the theoretical isochrone for x = 60 Myr, Z = 0.008, /j = 18.5, Eb-v = 0.2. Empirical 
colors for Cepheids (Madore & Freedman 1991) are shown with triangles and mark the location 
of the instability strip in the diagram. Stars show fiducial color sequence of supergiants (I-II) from 
W92. Dashed line shows the tip of the ZAMS at = 18.5, corresponding to hottest 03-06 dwarfs. 
Colors of Cepheids, supergiants and O dwarfs are unreddened. The reddening vector corresponds 
to Eb-v — 0.2. 
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Fig. 7. — Part of the CMD showing intermediate and old LMC populations. Theoretical isochrones 
represent prototype populations — solid line: x = 11 Gyr, Z = 0.004, fj = 18.5, Eb-v = 0-2; dashed 
line: x = 4 Gyr, Z = 0.004, /j = 18.5, Eb-v = 0-2. The long-dashed line shows the isochrone for 
Galactic RGB stars, x = 9 Gyr, Z = 0.019, /j = 10.0, Eb^v = 0.0. Fiducial unreddened RGB color 
sequence (W92) at /u = 18.5 is shown with squares, with approximate spectral types along the 
sequence marked in the top axis. The reddening vector corresponds to Eb-v = 0.2. 

the LMC sources in Region C have > 10.5. The similarity in the shapes of central isopleths 
between Fig. ^ and ^ suggests they are lower mass young supergiants with ages 300 — 500 Myr, 
evol ving into Region I (Figure |6|a). These stars trace the bar of the Cloud (G rebel & Brandner 
1998). Some contribution from more massive supergiants, including longer-period Cepheids 
(P < 100"^) may also be present. 



3.5. Region D: Disk G— M Dwarfs and LMC RGB and Early AGB Stars 

Region D is the most heavily populated area of the CMD: it includes more than a half of 
all sources in the sample. Because of its position in the CMD and the large color range it spans 
(0.25 < J — Ks < 1.2), this region is also the most inhomogeneous. The spatial distribution, 
shown in Figure 0D shows both the foreground and LMC populations (note the distorting effect 
of Galactic populations on outer LMC isopleths). The observed CMD is distinctly bimodal in this 
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region, with the red half populated by RGB and early AGB stars in the LMC, and the blue half 
populated mostly by G — M dwarfs in the Galaxy. The AGB stars in the red half of the region 
{J — Ks > 0.7) are in their 'early- AGB' (E-AGB) phase, during which the energy is produced in 
the thick helium shell and outer hydrogen shell is extinguished. These stars have recently passed 
the base of the AGB, the so-called 'AGB-bump' at Ks ^ 16 that marks the transition from core 
to shell helium burning ( |Castellani a/. 1991] ). The AGB-bump was first observed by Hardy et 
al. (1984) in their CMD of the LMC bar. At only ^ 1 mag brighter than the horizontal branch, 
this feature is not visible in the CMD in Figure |] although it is present in deeper data (Figure [TQ| ). 
Empirically, most stars at the E-AGB are M type (oxygen-rich). 

Only foreground stars are a minor contributor to the red half of Region D. (cf. Figure 
Galactic dwarfs in this region have ;U~8 — 11 (r~ 0.4— 1.6 kpc). Populations contributing 
insignificantly to the source density in this region include young supergiants, Cepheids, 
intermediate mass red stars in the vertically extended red clump (VRC; see § P.10| ). 



3.6. Region E: Upper RGB and Tip of the RGB 

Region E covers the upper RGB and includes the tip of the RGB (see §Q). Most of these stars 
are on the first-ascent red giant branch; they have degenerate helium cores and hydrogen burning 
shells. The majority of these stars have ages anywhere between 1 and 15 Gyr old. The tip of the 
RGB is defined the helium flash, the ignition of the degenerate helium core in old (low-mass) 
stars (Renzini & Fusi Pecci 1989). Stars at the TRGB ignite helium in their cores and evolve 
rapidly to the horizontal branch. The region also contains a significant fraction of AGB stars in 
transition from E-AGB to TP- AGB, the stage at which the outer hydrogen shell is re-ignited (I ben 
& Renzini 1983). During thermal pulses, the star begins alternating between hydrogen and helium 
shell burning. The transition from E-AGB to TP-AGB is theoretically predicted to occur near 
the TRGB. While on the TP-AGB, these stars may also experience the shorter-term atmospheric 
puls ations that lead to Mira-type variability. Analysis of MACHO data (|Alves et al. 1998^ W ood 



1999^" suggests that essentially all stars brighter and redder than the TRGB are variable. Most of 
the E-AGB stars in this region are M stars. Extrapolated to brighter magnitudes, the sequence of 
oxygen-rich AGB stars extends to Regions F and G (§§^]7|, |3^ . 

Stars in this region of the CMD carry the most weight in our analysis of the RGB-i-AGB 
luminosity function (see Their spatial distribution is relatively smooth, showing strong disk 
and bar components. Note the absence of significant foreground population in Figure the 
outer contours are elliptical in shape. A small fraction of foreground sources in this region is due 
to disk M dwarfs. Their density is steadily increasing toward fainter magnitudes (cf. Figure pb. 
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3.7. Region F: O-Rich AGBs 

Region F contains primarily oxygen-rich AGB stars of intermediate age ( > 1 Gyr) that are 
the descendants of stars in Region E (note the similarity between Fig. ^ and These are 
E-AGB and TP- AGB stars. The outer CMD isopleth in this region (Figure ^ is distorted and 
extends into Region J and indicates the presence of carbon stars in this region. During thermal 
pulses, the outer convective envelope may reach into the region where He has been transformed 
into C and bring carbon-enriched material to the surface. This dredge-up process leads to an 
increase in C/0 ratio, and M stars in Regions F and G may become carbon stars. In the CMD, 
carbon stars form a 'branch' with redder colors, J — Ks> lA (see § |3.1 1| ). Some fraction of 



Region F stars are LPVs (see Figure]^ and reddened supergiants. Figures ^ and G do not show 
isopleths due to a Galactic population. The Galactic component in these regions are negligible 
because it is both too bright and too red for disk M dwarfs and too faint for Galactic AGB stars. 



3.8. Region G: AGB Stars 

Region G contains the most massive stars with degenerate C/0 cores. This is a population of 
young AGB, post core-helium burning stars with initial masses between about 5-8 solar masses. 
These are too short-lived to become carbon stars (0.1 — 1 Gyr old), but not massive enough to 
become red supergiants. Similar to Region F, this region also includes LPVs (Figure |8|). The 
period-color relation for oxygen-rich Miras derived from Feast et al. (1989) rather closely traces 
the young AGB branch of the CMD outlined by Regions F and G ( [Weinberg & Nikolaev 2000[ ). 
Region G lies at bright enough apparent magnitudes that the foreground density of M dwarfs is 
low, even in comparison to the relatively small number of 0-rich luminous AGB stars in the LMC. 



3.9. Region H: LMC K— M Supergiants, Galactic M Dwarfs, K— M Giants 

Panel H of Figure ^ reveals an LMC population. The spatial distribution is similar to the 
distribution of young OB stars (Region A), suggesting that these objects are also relatively young. 
At/j= 18.5, these stars are too bright to be normal M giants. Based on their near infrared colors, 
we identify the Region H sources as supergiants of M type. They trace the spiral structure of the 
LMC and do not show significant overdensity in the bar of the Cloud, consistent with a young 
popu lation. The masses of these stars are believed to be ^ 2 — 9 solar masses ( Bertelli et al. 
1985). In the evolutionary sequence, these stars are descendants of stars in Region A (note the 
similarity between corresponding panels in Figure These stars are also the high mass extension 
of the VRC. 
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4 




Fig. 8. — Portion of the CMD showing evolved stars. Crosses indicate M supergiants Ml — M4 
from Elias et al. (1985), solid triangles — K-type LPVs, open squares — M-type LPVs, and solid 
squares — C-type LPVs from Hughes & Wood (1990). Arrows at left show the theoretical tip of 
the RGB stars and the lower luminosity limit for thermally pulsing AGBs. The reddening vector is 
fox Eb-v = 0.5. 

In Figure ^ we plot the observed colors of Ml — M4 supergiants from the sample of Elias et 
al. (1985). The colors of their sample occupy a portion of Region H, supporting our identification. 
The Galactic foreground consists of roughly equal contributions from disk K — M giants and M 
dwarfs, but their overall contribution to the source density in the region is only a few percent. This 
is confirmed by the absence of the Galactic isopleths in Figure ^i. 
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3.10. Region I: LMC Intermediate-Mass Red Supergiants, Galactic K — M Dwarfs 

Region I is located at the center of the CMD, at 0.7 <J — Ks< 1 .0. The observed overdensity 
in the CMD is associated with the vertically extended red clump (Zaritsky & Lin 1997). This 
feature consists of intermediate mass stars and is the low-mass extension of the red supergiants 
(Region H). The VRC extends upward from the red clump at ~ 17 and becomes visible in the 
CMD near Ks — 13.5. At the this point, the redward slope of the RGB is sufficient to distinguish 
the VRC. 

The spatial distribution is dominated by the bar and shows traces of the spiral structure. 
We conclude that this LMC population is young, with the age < 500 Myr. The major LMC 
contributors to the source density in this region are K and M supergiants. This is supported by 
the overall similarity of LMC isopleths in Figures ^ and and also the fact that Region I is at 
the extension of Region H to fainter magnitudes and lower masses. Figure |^ shows K and M type 
Miras and SR variables in the sample of Hughes & Wood (1990). A significant fraction of their 
variables falls in this CMD region suggesting that some of these 2MASS stars also are variables. 

The distribution in Figure |^I also reveals foreground populations. The Galactic foreground 
consists of M and late K dwarfs. Galactic giants contribute less than 5% of the foreground density. 
The dwarfs are located in the disk, with distance moduli /j = 5 — 8(r~0.1— 0.4 kpc; Figure 
The contribution from the Milky Way halo is smaller than 0.0005% by number for this and for all 
other regions of the CMD. 



3.11. Region J: Carbon Stars in the LMC 

At J — Ks > 1.4, Region J sources are primarily carbon-rich TP-AGB stars. These stars are 
descendants of oxygen-rich TP-AGBs in Regions F and G. Their outer layers are enriched in C 
through convection from stellar interior. As mentioned in § p.6| , most of these stars are long-period 
variables. The variability cannot be determined based on single epoch 2MASS data, but the 
well-defined sequence motivates a follow-up campaign. Figure |^ shows the sample of C-rich 
LPVs from Hughes & Wood (1990) overplotted on the 2MASS CMD. The contamination by 
M-type LPVs is small. The spatial distribution of C stars in the field is similar to the distribution 
of their precursors (Figs. ^). The distribution is rather smooth and shows a loop of stellar 
materiafl, which has been described by Westerlund (1964). The loop is the extension of the main 
northern spiral arm circling the main body of the system and returning toward the bar after a 
nearly complete turn. 



^The feature to the SE of the bar in Figure^, near a « 5''50"', 5 « —73° represents a hole in the disk. 
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Sources in this region of the CMD offer the best opportunity to study the three-dimensional 
structure of the LMC for two reasons. First, the spatial coverage of the Cloud achieved by 2MASS 
is total and allows to probe the entire LMC. Second, as long-period variables. Region J stars are 
potentially good standard candles, since their intrinsic luminosity can be characterized based on 
their period or color. Given the selection efficiency^ and easily quantifiable intrinsic brightness 
through the period-luminosity-color relation (e.g.. Feast et al. 1989), these stars are excellent 
probes of the LMC structure along the line of sight. Preliminary results (WN) indicate that the 
width of the intrinsic brightness distribution is smaller than Om = 0.2 magnitudes in a narrow 
color range, A(7 — ~ 0. 1 . At this accuracy, these standard candles can resolve features in the 
LMC at Ar ~ 4.5 kpc. 2MASS detected approximately 10"^ potential carbon LPVs and these are 
sufficient to attain a reasonable confidence level in the inferred spatial structure. In Weinberg & 
Nikolaev (2000), we present our study of the three-dimensional structure of the LMC. 



3.12. Region K: Dusty AGBs 

An extension of Region J, Region K contains extremely red objects. We identify them 
with obscured AGB carbon-rich stars. Their large J — Ks colors are due to dusty circumstellar 
envelopes (Eb-v ^1)- The latter is confirmed by the appearance of their spatial and CMD 
distributions: (1) Figure ^ shows traces of the spiral structure outlined by these sources; and 
(2) the distribution in the CMD spreads from the end of Region J in the direction of reddening 
vector. Matching with existing near infrared photometry of obscured AGB stars in the LMC 
( [Zijlstra et al. 199^ ; |van Loon et al. 199B| ) shows that most of these sources are indeed in this 



region of the CMD. Other extremely red populations could also be found here, e.g. 'cocoon' stars 
( IReid 1991D , or OH/IR stars ( |Wood a/. 1992| ; |van Loon aZ. 1998D . In addition, two of the 



known LMC protostars, N159-P1 and N159-P2 ([Tones et al. 19860 also fall in this region. 



3.13. Region L: Reddened LMC M Giants, Galactic M Dwarfs and 2MASS Galaxies 

Stellar sources in Region L are reddened M giants in the LMC and a small number of 
reddened Galactic M dwarfs. However, a significant number of sources in this region are 
background galaxies. The predicted CMD density in Region L due to Galactic stars is too low 
even after the decrease in the photometric quality near the flux limit has been taken into account; 
the decreasing signal-to-noise ratio causes the apparent widening of the contour levels (see 



owing to their extremely red colors, these stars are uncontaminated by other populations 
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Figure According to Jarrett (1998), more than 90% of 2MASS galaxies have colors redder 
than 7 -if., = 1. 

The spatial distribution of sources (Figure ^) shows some overdensity toward the center of 
the LMC. The densest part of the diagram corresponds to the position of 30 Doradus complex. 
Traces of spiral structure of the LMC are also visible. Based on their colors, these sources are 
heavily obscured RGB stars in the LMC: they lie in the direction of the reddening vector from the 
RGB. The inferred reddening for these sources, Eg^y ~ 0.5, is consistent with the extended tail 
of the LMC reddening distribution ( Parris et al. 1997| ). Region L also includes contribution from 



massive (> IOM0) protostars and ultra-compact H II regions, see Figure |S|. 

A population of dwarfs is implied by the outer isopleths of Figure ^ that show the increase 
in the direction of Galactic center. These are local M dwarfs in the disk of the Milky Way, with 
Ai~5-8 (r~ 100-400 pc). 



4. Luminosity Function of LMC RGB and AGB Populations 

We derive the LMC giant branch luminosity function (LF) from the color-magnitude diagram 
after subtracting Galactic foreground. Since we did not have the access to Galactic data at the 
time, the foreground contribution was estimated from three small fields located at the edges of 
our LMC field. We then scaled the resulting foreground CMD to the entire LMC field by using 
the estimate for the number of Galactic sources from our synthetic model. Figure ^ shows 
the field CMD for LMC populations only, after subtracting Galactic foreground. The expected 
Galactic source counts is ~ 4 x 10^ or about 50%. The uncertainties in the observed CMD and in 
the Galactic model result in negative density regions (dotted contours) in Figure |^. The average 
negative density is — lO-^'^ mag^^. 

The luminosity function of the LMC giants is obtained by projecting the color-magnitude 
diagram perpendicular to the giant branch ridge line. The function is normalized to unity. 
Numerical values for the RGB luminosity function are given in Table |l] and shown in the inset 
to Figure ^ A strong feature of the LF is a significant excess at Ks ~ 12.5, due to the TRGB. 
From the analysis of the derivative of the apparent luminosity function, we derive the position 
of the TRGB at iS:^ = 12.3 ± 0. 1 . Brightward of the TRGB, the number of RGB stars drops off. 
The increase in the number density at the faint end, Ks> 14, is due to the increased contribution 
from Galactic M dwarfs (cf. Fig.^. At the bright end of the magnitude range, 11 < < 12, 
the luminosity function is nearly constant. It is also well above the expected number from 
extrapolated RGB counts. The fraction of RGB stars is small in this magnitude range and most of 
the stars contributing to the LF are on the AGB. As discussed above (§pJ[), these stars tend to be 
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oxygen-rich, but carbon-rich AGBs (and LPVs) are also present. 

We select two 2° x 1° fields, one near the optical center of the bar, and the other one 
at a = 93°, 6 = —67.5° (J2000.0), near the outer loop, to compare the observed M giants 
luminosity functions in distinct LMC environments. The outer field is probing the LMC's 



outer loop delineated by the evolved stars (see § |3.1 1| ). For each of the fields, we subtract the 
estimated Galactic foreground density scaled to the surface area of the fields. Similar to Fig. 
inaccuracies in both observations and models produce negative density regions. However, the 
average negative density in these regions is small compared to the giant branch, only — lO*^ "^ 
mag^^ for bar region and — lO'^'^ mag^^ for loop field. Comparing the two CMDs qualitatively, 
we note that contribution of young OB stars and supergiants, at J — Ks < 0.2, appears stronger 
in the central regions of the LMC. Even though the luminosity functions in the respective fields 
appear different, careful analysis shows that the difference is superficial: %^/d.o.f.^ 0.2, which 
gives no motivation to entertain any difference in the parent populations. The bar LF is similar to 
the luminosity function of the entire field (Figure ^), which is not surprising because the source 
density in the LMC field is dominated by the LMC bar. The bar luminosity function also has a 
pronounced TRGB at ~ 12.3 and excess number density due to AGBs at 11 < < 12. The 
sharp increase due to Galactic M stars, seen in Figure ^ at Ks > 14, has disappeared in the bar 
field, because we have boosted the ratio of LMC/Galactic counts by narrowing down the field to 
the area of greatest LMC density. The off -bar LF shows only a mild increase in the source counts 
at the location of TRGB, but has the same, roughly constant profile at Ks < 12, due to the AGB 
population, visible in the other two luminosity functions. To quantify both LFs, we present their 
numerical values in Table [I]. The luminosity functions are given in relative units, normalized 
to unity. The table also gives the source counts for the LMC giant branch. For the entire LMC 
field we present the total counts per magnitude bin, and for the two smaller fields we give stellar 
density (counts mag^^ deg^^). 



We fit theoretical isochrones ( pirardi et al. 2000| ) to each giant branch in Figure ^ to test for 



differences in metallicity between central and outer parts of the LMC. We chose 20 equally-spaced 
grid points in Ks magnitude between 14.3 and 12.3 and compute the peak in the distribution in 
J — Ks at these fixed Ks points. The difference between an isochrone (model) and the RGB (data) 
is characterized by the cost function (mean integrated square error): 

/ = Y}^J~Ks)j,RGB - {J-Ks)j,,so?+w[Kj''''^ - 12.3]2, (1) 
j 

where the second term is weighted measure of the match between theoretical magnitude at helium 
flash and observed TRGB. The weight w is an adjustable parameter on the order of unity. The 
cost function (|l]) is minimized on a grid of parameter values, where the free parameters are the 
log-age X, metallicity Z, distance modulus /j and average reddening Eg-v- The best fit isochrones 
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Fig. 9. — Color-magnitude diagrams of the LMC populations and the apparent luminosity 
functions of the LMC giants. Panels show CMDs of (a) entire field, (b) bar field, and (c) outer loop 
field (see text). The CMDs are constructed by subtracting Galactic foreground contribution from 
CMD in Figure |3|, normalized to the same sky area. Contour levels in all panels are logarithmic, 
spaced by 0.5, from 2.5 to 5.5 (a), 2.0 to 4.0 (b), and from 1.0 to 3.5 (c). Dotted lines indicate 
negative density regions. Luminosity functions, normalized to unity, are shown in the insets of 
each panel. 
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Table 1 : Apparent luminosity function and source number density for the LMC giants 



Ks 


4> (mag 




log number density 










(mag 1 


deg-2) 






LMC 


Bar 


Loop 


LMC'' 


Bar 


Loop 


10.8 


0.03 


0.04 


0.05 


3.16 


2.17 


1.16 


11.0 


0.05 


0.06 


0.09 


3.44 


2.32 


1.39 


11.2 


0.07 


0.08 


0.10 


3.61 


2.50 


1.46 


11.4 


0.08 


0.09 


0.09 


3.65 


2.54 


1.40 


11.6 


0.08 


0.09 


0.08 


3.66 


2.54 


1.35 


11.8 


0.08 


0.08 


0.08 


3.64 


2.48 


1.33 


12.0 


0.09 


0.12 


0.11 


3.68 


2.63 


1.47 


12.2 


0.19 


0.24 


0.21 


4.04 


2.96 


1.76 


12.4 


0.29 


0.31 


0.26 


4.22 


3.06 


1.87 


12.6 


0.36 


0.37 


0.30 


4.31 


3.14 


1.92 


12.8 


0.40 


0.40 


0.35 


4.35 


3.17 


1.99 


13.0 


0.40 


0.41 


0.36 


4.36 


3.18 


2.00 


13.2 


0.43 


0.42 


0.36 


4.39 


3.20 


2.01 


13.4 


0.44 


0.43 


0.45 


4.39 


3.20 


2.11 


13.6 


0.44 


0.45 


0.47 


4.40 


3.23 


2.14 


13.8 


0.47 


0.47 


0.46 


4.43 


3.25 


2.13 


14.0 


0.54 


0.49 


0.57 


4.48 


3.26 


2.22 


14.2 


0.58 


0.47 


0.61 


4.51 


3.24 


2.25 



'^The luminosity functions are given in relative counts per magnitude. 
''Number density for the entire field is given as source counts per 
magnitude bin. 



are as follows: {x,Z,iu,Eb-v) = (9.8 ±0.3,0.004ljJ;[j[;?, 18.45 ±0.11,0.21 ±0.07) for the central 
field and {x,Z,^,Eb-v) = (9.8 ±0.4, 0.004+5];gj;2, 18. 50±0.13, 0.13 ±0.09) for the outer field 
(all errors statistical). These results imply an age range for RGB populations from 3 to 13 Gyr 
with an average of 6 Gyr. The slope degeneracy of the isochrones in the RGB makes specific tests 
of star formation history difficult. In particular, our preliminary RGB isochrone analysis cannot 
distinguish between a continuous and single/multiple burst star formation history of the Cloud 
prior to approximately 4 Gyr ago. Overall, our results do not indicate a radial metallicity gradient 
and provide only marginal evidence for larger reddening in central fields. The absence of strong 
metallicity gradients in the LMC is in agreement with results of Olszewski et al. (1991), who 
found no evidence for abundance gradient for cluster system. Constant C/M star ratio across the 
face of the LMC ( [Westerlund 1997[ ) and Cepheid abundances ( parris et al. 1983| ) also support 
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this result. Our results imply the range of abundances for field populations —0.8 <[Fe/H]< —0.5, 
which is in good agreement with the mean abundance [Fe/H]= —0.58 ±0.05 (systematic) ±0.30 
(statistical) for the inner LMC disk ( |Cole 1999| ). Our results agree with the disk abundance 
[Fe/H]= -0.7 dCowley & Hart wick \ W^ , and with resuks of Bica et al. (1998), who derived the 



range —1.1 <[Fe/H]< —0.4 with the average < [Fe/H] >= —0.61 ±0.11 from fields in the outer 
disk of the LMC. 

The adopted detection threshold (SNR = 10; §0) leads to the effective completeness limit 
of Ks ~ 14.3 in our data (cf. Figure |]). With this flux level, only the upper RGB is visible, 
leaving out both the AGB-bump at ~ 16 and the red clump dX K^'^ 11 . To resolve the giant 
branch down to 16 — 17, we use 2MASS engineering data, which includes six LMC scans, 
positioned as shown in Figure [T^. Each of the 'deep' scans has six times the standard exposure. 
The color-magnitude diagram of the deep data is shown in Figure [TT|. Total number of sources in 




BHO" fii-OO" 5''20" 4H0" 



R.A. (J2000.0) 



Fig. 10. — Positions of six 'deep' scans in the sky. Each scan has six times the standard integration 
time. 



the diagram is 87,093, of which 69, 878 (~ 80%) are in the three bar scans. Because of the bar 
dominance in deep data, the CMDs of the entire deep sample and the bar scans (Figs. [TT] a and 
pT|b, respectively) are similar. The increased sensitivity reaches the AGB-bump at 7 — = 0.7, 
Ks — 15.8, but still shy of the red giant clump. As with the main dataset, we quantify the deep 
RGB population by fitting isochrones. The resulting best fit parameters are Z = 0.004^q qqj, 
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Table 2: Regions of interest in the 2MASS CMP of the Large Magellanic Cloud. 



Region 


Nsrc 


fCaf 


Region boundaries 


Dominant Spectral Types'^ 


A 


6,659 


0.15 


l\<Ks< \4J5,J-Ks<0.2 


B-A I II, 03-06 V 


B 


77,204 


0.80 


5.5<K,< 13.5,0.2 < J -K,< 0.5 


F-KV 


C 


62,713 


0.80 


5<K,< 13.5, 0.5 < < 0.8 


K V, K III 


D 


440,472 


0.45 


13.5 <K,< 14.75, 0.2 <7-i^,< 1.2 


K-M III, F-M V 


E 


166,263 


0.05 


12<i^,< 13.5,0.9<7-i^,< 1.2 


M III, M V 


F 


22,134 





10.5<i^,<12, l<7-i^,<1.3 


M, MS 


G 


1,438 





8<i^,< 10.5, 1.2<7-i^,< 1.5 


M, MS 


H 


2,450 


0.05 


7<is:,<ll, l<7-i^,<1.3 


MI-II 


I 


21,986 


0.55 


ll<if,< 13,0.75<7-i^,< 1 


K-M I-II, K-M V, M III 


J 


8,229 





9.5<i^, < 11.5, 1.4<7-i^,<2 


cm 


K 


2,212 





9<if,< 13,2<7-iC,<5 


cm 


L 


8,940 


0.01 


12.5 <K,< 14.75, 1.2 < 7-i^, < 2.5 


M late V 



''Fraction of Galactic sources estimated from synthetic W92 model 
''Based on 7 — Ks color and W92; LMC populations in boldface 



X = 9.7 ±0.3, fi= 18.50 ±0.12, Eb-v = 0.19 ±0.08. The uncertainties here are statistical errors, 
derived from the shape of the %^ surface near minimum. The range of ages for RGB populations 
inferred from the deep data is similar to that derived for the main data set: from 3 to 10 Gyr, with 
the average of 5 Gyr. The other parameters are also consistent with the values derived from the 
regular 2MASS data. These estimates are in good agreement with recent results in the literature, 
e.g., the average reddening Eg-v 
A' 



0.20 dHarris et al. 1997| ); the LMC distance from Key Project, 



18.5 (Mould et al. 1999D; the LMC distance from TRGB, ^ 



1999). 



18.59 ±0.09 ( Sakai et al. 
I 



5. Summary 

In this paper, we analyzed the near-infrared CMD of the Large Magellanic Cloud and 
identified the major stellar populations. The populations are identified based on isochrone fitting 
and matching the theoretical CMD colors of known populations to the observed CMD source 
density. Tables ^ and |^ summarize the contents of the CMD regions. The main points of this 
preliminary analysis of 2MASS data are the following: 

• The quantity and the quality of 2MASS data allow unprecedented look at the entire 
LMC. 2MASS has produced a rich sample of LMC sources, a few million stars, with the 
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Fig. 11. — 'Deep' color-magnitude diagrams of the LMC. Panels show CMDs of (a) all six deep 
scans combined, (b) three bar scans, and (c) three outer field scans. Panel (d) is the same as panel 
(a), except also shows the best-fit isochrone for the RGB. Galactic foreground is not subtracted. 
The contour levels are spaced logarithmically by 0.5, from —3.6 to —0.1 (a,b,d), and from —3.7 
to —0.2 (c). The lower RGB is enhanced as compared to Figure 9. Weakness of RGB in panel (c) 
indicates strong relative contribution from Galactic foreground. The red giant clump is just below 
the completeness limit in this diagram, at ~ 11, J — ~ 0.65. 
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photometric accuracy of 3-4%. JHK 2MASS photometry is potentially useful for studying 
the star-formation history of the Cloud. Cross-correlating 2MASS database with existing 
catalogs will provide homogeneous and accurate IR photometry of supergiants (Sanduleak 
catalog), Wolf-Rayet stars (Breysacher catalog), Cepheids (OGLE and EROS datasets), 
LBVs and LPVs; 

• The color-color diagram is generally ill-suited to distinguish between giant (HI) and dwarf 
(V) populations, especially in the color range 0.5 < 7 — < 0.8. Nevertheless, the diagram 
may be useful in identifying some candidate LMC objects with infrared excess, such as 
obscured AGE stars, B[e] stars, or LMC protostars. In addition, the distribution of 7 — 
colors of a population in a narrow J — H color range is a sensitive reddening test; 

• Major populations can be identified based on the comparison of observed CMD features 
with theoretical positions of known populations. Isochrone overplotting provides tentative 
age and metallicity estimates. We identify a substantial LMC population of AGBs ( > 10^ 
sources), and obscured AGBs (~ 2000 sources); 

• The luminosity function of the LMC giants is determined and tabulated. We find the RGB 
tip at = 12.3 ± 0. 1 . Our preliminary analysis of luminosity functions in two test fields 
suggest that luminosity function is the same in the bar and the outer regions of the Cloud; 

• Fitting isochrones to the location of the giant branch (including TRGB) gives metal 
abundances consistent between fields. In particular, we derive average metallicity 

Z = 0.004;^q"qqj for our fields. Analysis of deep data gives the same average metallicity. 
Our results confirm the absence of strong radial metallicity gradient in the field populations 
of the LMC; 

• The estimates of the distance modulus obtained from our isochrone fits to the RGB, 
are consistent with each other and the most recent results in the literature. The average 
reddening is marginally different between the bar and the outer field, the Eb-v for the bar 
field being greater. Distance modulus and reddening estimates from the analysis of deep 
data produces similar values; 

• The ages of dominant RGB populations fall in the range from 3 to 13 Gyr, with the average 
age ~ 6 Gyr. Isochrone fits to the deep data produce similar age interval, from 3 to 10 
Gyr, with the average age ~ 5 Gyr. A more detailed isochrone analysis is required to draw 
conclusions about the history of star formation in the LMC prior to 3 — 4 Gyr ago; 

• Carbon-rich long-period variables are noted as potential standard candles. Due to their 
significant numbers and narrow luminosity range (which may be parametrized through a 
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period-luminosity or luminosity-color relations), these stars are ideally suited for studying 
the structure of the LMC along the line of sight. 
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Table 3: LMC stellar populations in the 2MASS CMD 



Stellar types Regions^ Typical age 

Very Young 

Centrally concentrated, localized to star-forming regions, trace spiral structure, weakly trace bar 

03-06 dwarfs A < 10 Myr 

Red supergiants, 5 - SMq H < 50 Myr 

Luminous AGB stars, 0-rich LPVs, M^5- SM© G, H, F ~ 40 - 100 Myr 

Blue and yellow supergiants, LBVs A, B, C < 100 Myr 

Massive protostars or cocooned OB associations L < 5 Myr 

Young 
Bar-dominated 

Luminous E-AGB stars G, H 200 - 800 Myr 

Core He-burning giants and supergiants, M ~ 2 - 5Mq, LPVs I, D, C, H 100 - 900 Myr 

E-AGB, oxygen-rich LPVs F, E, D < 1 Gyr 

Carbon stars, C-rich LPVs J, K < 1 Gyr 

Intermediate and Old 

Disk and bar 

Low- and intermediate-mass RGB stars D, E, L 1 — 15 Gyr 

0-rich AGB stars, M-S-C stars F 1 4 Gyr 

C-rich TP- AGB stars J 1 - 4 Gyr 

Dust-enshrouded TP- AGB, carbon stars K 1-4 Gyr 

Foreground and Background 
Disk populations and extragalactic component 

Disk main-sequence tumoff stars B < 7 — 9 Gyr 

Nearby K dwarfs, red clump and red HB stars C, I < 9 Gyr 

Local F-M dwarfs D, E varies 

Background galaxies L 



''Where a population is the dominant contributor to a region, it is labeled in bold type 



